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Abstract 

We study a simple microscopic model for the one-dimensional stochastic motion 
of a (non) relativistic Brownian particle, embedded into a heat bath consisting of 
(non) relativistic particles. The stationary momentum distributions are identified 
self-consistently (for both Brownian and heat bath particles) by means of two cou- 
pled integral criteria. The latter follow directly from the kinematic conservation laws 
for the microscopic collision processes, provided one additionally assumes proba- 
bilistic independence of the initial momenta. It is shown that, in the nonrelativistic 
case, the integral criteria do correctly identify the Maxwellian momentum distribu- 
tions as stationary (invariant) solutions. Subsequently, we apply the same criteria 
to the relativistic case. Surprisingly, we find here that the stationary momentum 
distributions differ slightly from the standard Jiittner distribution by an additional 
prefactor proportional to the inverse relativistic kinetic energy. 
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1 Introduction 



The implementation of the Brownian motion concept [1-6] into special relativ- 
ity [7,8] represents a longstanding issue in mathematical and statistical physics 
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(classical references are [9-11]; more recent contribntions include [12-22]; for 
a kinetic theory approach, see [23-25]). In two recent papers [26,27] we have 
discussed in detail how one can construct Langevin equations for relativistic 
Brownian motions (see Debbasch et al. [28, 29] and Zygadlo [30] for similar 
approaches, and also Dunkel and Hanggi [31]). Thereby, it was demonstrated 
that, in general, the relativistic Langevin equation per se cannot uniquely de- 
termine the corresponding Fokker-Planck equation (FPE). This dilemma is 
caused by the fact that the relativistic Langevin equations, e.g., if written 
in laboratory coordinates, may exhibit a multiplicative couphng between a 
function of the momentum coordinate and a Gaussian white noise process 
(laboratory frame = rest frame of the heat bath). Thus, depending on the 
choice of the discretization rule, one obtains different forms of relativistic 
FPEs characterized by different stationary solutions. 

In Refs. [26, 27] we have analyzed the three most popular discretization rules 
for Langevin equations with multiplicative noise, namely, Ito's pre-point dis- 
cretization rule [32,33], the Fisk-Stratonovich mid-point rule [34-37], and 
the Hanggi- Klimontovich (HK) post-point rule [38-41]. As a main result it 
was found that only the HK interpretation of the Langevin equation yields a 
FPE, whose stationary solution coincides with the one-dimensional relativistic 
Jiittner-Maxwell distribution, as known from Jiittner's early work on the rel- 
ativistic gas [42,43] and also from the relativistic kinetic theory [23,24]. Thus, 
in absence of other qualifying criteria, one may conlude that the post-point 
discretization rule is favorable. However, it naturally arises the question, if one 
can gain additional insights by studying the microscopic collision processes, 
which cause the stochastic motion of a Brownian particle. 

The present paper intends to partially address this question for the one- 
dimensional (ID) case; i.e., instead of focusing on the 'macroscopic' Langevin 
formulation, we will study here Brownian motions by means of a simple mi- 
croscopic ID model. Since, on the level of the Brownian motion approach, 
the exact nature of the microscopic interaction forces is usually negligible we 
shall restrict ourselves to purely elastic collisions between a Brownian parti- 
cle and the constituents of a surrounding heat bath. As will be demonstrated 
below, this simple model suffices to identify stationary (invariant) momentum 
distributions for both non-relativistic and relativistic ID collision processes ^ . 

The paper is organized as follows: Section 2 is dedicated to the non-relativistic 
case, serving as the test example for our approach. After briefly summarizing 
the basic model assumptions, we will derive two coupled integral criteria (Sec. 
2.3), which can be used to identify the invariant momentum distributions for 



More precisely, we shall additionally require the probabilistic independence of the 
initial momenta during each single collision - but this is a rather reasonable, weak 
restriction. 
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the heat bath particles and the Brownian particles, respectively. It is shown 
that these integral criteria do indeed yield the correct stationary momentum 
distribution of the non-relativistic Brownian motion, namely, the well-known 
MaxweUian momentum distribution (Sec. 2.4). Afterwards, the same method 
is applied to the relativistic case (Sec. 3). Remarkably, we find that, under 
exactly same preconditions, the invariant relativistic distributions do not ex- 
actly correspond to the standard Jiittner distributions, but rather to modified 
Jiittner functions including an additional prefactor 1/E, where E is the rela- 
tivistic kinetic energy of the particle under consideration. The paper concludes 
with a summary and a discussion of the results (Sec. 4). 



2 Non-Relativistic Brownian motions 

In this part, we review non-relativistic ID-Brownian motions. Later on, we 
will pursue an analogous approach to identify the stationary momentum dis- 
tribution of a relativistic Brownian particle, embedded into a relativistic heat 
bath. 

2.1 Basic model assumptions 

In the three-dimensional (3D) case a simple idealized model for Brownian mo- 
tions can be imagined as follows: Consider a (infinitely heavy) box of volume 

V, possessing diathermal [44] walls and being at rest in the inertial laboratory 
frame Sq. Let this box contain a homogeneous, quasi-ideal (weakly interact- 
ing) gas, consisting of approximately point-like particles with identical masses 
m. Further, assume that the gas (or liquid) particles - referred to as 'heat 
bath', hereafter - surround a Brownian particle of mass M ^ m. Then, due 
to frequent elastic collisions with heat bath particles, the Brownian particle 
performs 3D random motions. Given the distribution of the heat bath, the 
stochastic dynamics of the Brownian particle is determined by the collision 
kinematics (cross-sections) governing the interaction with the heat bath par- 
ticles. 

If, as in the present paper, one wishes to study the ID case, slight modifi- 
cations of the above model are necessary. The reason is that, typically, two 
particles cannot simply exchange positions if their motions are confined to 
one dimension. To circumvent this problem, we shall therefore imagine the 
heat bath particles as having fixed positions on a ID lattice but non-vanishing 
momenta. Correspondingly, in this ID (lattice) model the Brownian particle 
may jump from one lattice point to the next during one time step. Addi- 
tionally, we will impose that at each lattice point there does indeed occur an 
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elastic interaction in accordance with the laws of momentum and energy con- 
servation. Mathematically, the latter assumption corresponds to considering 
distributions conditional on the event "a collision has occurred" . 

Based on this idealized ID (lattice) model, it is our primary objective to deter- 
mine self-consistently the invariant (i.e., stationary) momentum distributions 
for both heat bath and Brownian particles. To this end, we shall next summa- 
rize the simple kinematic equations governing the collisions in this model. By 
interpreting the momentum coordinates as coupled random variables, we then 
derive in Sec. 2.3 two general integral criteria which have to be satisfied by 
the stationary distributions. As shown in Sec. 2.4, in the non-relativistic case 
the stationary solutions are given by the well-known Maxwell distributions. 
The integral criteria apply to both non-relativistic and relativistic collisions; 
hence, we can use them later to also identify the stationary distributions of 
the corresponding relativistic model (cf. Sec. 3). 



2.2 Kinematics oj single collision events 



The momentum and energy balance per (elastic) collision reads 



E + e = E + e, (la) 
p + p=p + p, (lb) 

Here and below, capital letters refer to the Brownian particle and small letters 
to heat bath particles; quantities without (with) hat-symbols refer to the state 
before (after) the collision. In the non-relativistic case, we have, e.g., before 
the collision 



P = MV, p = mv, (2a) 

, e = ^, (2b) 

2M' 2m' ^ ' 

where v and V denote the velocities with respect to the laboratory frame Eq. 
Taking into account both conservation of momentum and (kinetic) energy, one 
finds for a single collision the elementary results 



P{p,P) 
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2m 

M + m 
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{M 
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\M + m 
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(3a) 
(3b) 
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We again stress that Eqs. (3) do implicitly assume that a collision has indeed 
occurred (otherwise, the momenta would remain unchanged); i.e., in mathe- 
matical terms, any results obtained by employing Eqs. (3) are valid conditional 
on the information that a collision event has taken place. 

Now let us suppose we know the joint two-particle PDF ip2ip, P) for the par- 
ticle momenta before the collision. Then, the kinematic laws (3) determine 
uniquely the marginal momentum PDFs $(-P) and after the collision, 
formally defined by 

$(P) = I dpMp,P), (4a) 
0(p)=|dPV'2(p,P), (4b) 

where ip2{p, P) is the joint momentum PDF after the collision (here and be- 
low integrals with unspecified boundaries range from — oo to -|-oo.) Our main 
objective will be to identify stationary (invariant) momentum distributions, 
satisfying hy definition 



Mp,p) = Mp,p)- (5) 

The latter condition just means that a stationary PDF must remain invariant 
in microscopic collisions. In particular look, we shall for particular stationary 
solutions which can be written in the product form 



V^2(p,P) = 0(p)$(P). (6) 

Mathematically, this corresponds to the assumption that, in the stationary 
state, the momenta P and p can be considered as independently distributed 
random variables. For stationary PDFs of the form (6) the stationarity crite- 
rion (7) reduces to 



$(P) = $(P), 0(p) = 0(p), (7) 

where ip2{p, P) = i'ip) ^{P) is the joint distribution after the collision. From 
a physical point of view, the independence assumption for (p, P) or, alter- 
natively, of (p, P) is guided by the experience that the well-known equilib- 
rium momentum distributions of quasi-ideal non-relativistic and relativistic 
N particle gases (i.e. the Maxwell and Jiittner distributions) can be written 
as products of one-particle momentum distributions. 

In order to be able to determine the stationary PDFs for a given collision 
kinematics, we next derive general integral criteria. It comes as no surprise 
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that, in the non-relativistic case, the stationary solutions will be given by a 
pair of Maxwellians. 



2.3 Integral criteria for stationary momentum distributions 



Consider two independently distributed random variables Y and Z, and a 
derived random variable X — X(Y,Z). The corresponding PDFs are denoted 
by ^x{X), ^viy) s^iid ^z{Z). The average of some test function g{X) with 
respect to $x can then be written as 



j dX g{X) ^x{X) ^ JdY jdZ g{X{Y, Z)) $y(y) ^z{Z). (8) 

Assuming that the (partially) inverse transformation Z = Z(Y,X) is well- 
defined (i.e., strictly monotonous for each fixed Y), and that Fubini's theo- 
rem [45] is applicable, we can rewrite the last equation in the form 



I dX g{X) ^x{X) ^jdY jdX 



dZ 



dX 

^jdX g{X) J dY 



g{X) ^y{Y)^z{Z{Y,X)) 
dZ 



dX 



^y{Y)^z{Z{Y,X)). (9) 



Since the latter equation holds for any test function g, one obtains the well- 
known transformation law 



dY 



dZ 



dX 



$y(y) ^z{Z{Y,X)). 



(10) 



For completeness, we note that Eq. (10) can also be obtained by starting from 



^x{x) - JdY I dZS{x - X{Y, Z)) ^y{Y) ^z{Z) (11) 

and performing the Z-integration (with 6 denoting the Dirac delta-function). 
We next consider an explicit example, which will be investigated in detail in 
the remainder. 



Example: P = P{p,P) and p = p{P,p). The idea is that we express the 
final momenta in terms of the initial momenta, cf. Eq. (3). That is, setting 
X = P,Y = p, Z = P and, correspondingly, $x = $, $y = 0, $z = $, we 
can write Eq. (10) as 



6 



HP) 



dp 



dP 



dP 



<i>{p)<^{P{p,P)). 



(12a) 



For a given pair of initial distributions ($,0), this equation can be used to 
calculate the PDF of the Brownian particle after the collision. Analogously, 
by setting X^p,Y^P,Z^p and $x = 0, $y = = 0, we find for 

the PDF of the heat bath particles 



0(p) 



dP 



dp 
dp 



$(P) 0(p(P,p)). 



(12b) 



In particular, for stationary distributions satisfying Eqs. (7) wc have 
and $ = $ and, hence, obtain from Eqs. (12) the integral criteria 



r dP 
^P)^J dp 



dP 
dp 
dp 



^{p) ^P{p,P)), 
$(P) 0(p(P,p)). 



(13a) 
(13b) 



Given a certain microscopic kinematic law, any pair ($, 0) satisfying Eqs. (13) 
provides a self-consistent stationary distribution. For completeness we men- 
tion that mathematically equivalent criteria are obtained by exchanging the 
positions of p and P as functional arguments, i.e., by considering P = P{P,p) 
and/or p = p{p, P), respectively. 

Before discussing solutions of Eqs. (13) for the non-relativistic Brownian mo- 
tion, it is worthwhile to stress the following fact: Since the derivation of 
Eqs. (13) is based on rather general assumptions, these integral equations 
can be applied to find the stationary PDF $ not only in the non-relativistic 
but also in the relativistic case (as will be done in Sec. 3). The additional 
mathematical assumption, underlying the derivation of Eqs. (13), is that the 
initial momenta P and p can be viewed as independently distributed random 
variables; i.e, loosely speaking, this postulate is the only point leaving some 
freedom for potential modifications, all other parts are dictated by the phys- 
ical conservation laws. As stated before, from a physical point of view, the 
independence assumption for (p, P) or, alternatively, of (p, P) is guided by 
the experience that the well-known equilibrium momentum distributions of 
quasi-ideal non-relativistic and relativistic particle gases (i.e. the Maxwell 
and Jiittner distributions) can usually be written as products of one-particle 
momentum distributions. 
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2-4 Stationarity of the Maxwell distribution 



In the last part of this section, we briefly outline that the integral criteria 
Eqs. (13) are satisfied by the normahzed Maxwell distributions 



$(P) = ( 



Hp) 



1/2 



exp 



2nmkBT 



2MkBTJ' 

1/2 / „2 

exp I 



2mkBT r 



(14a) 
(14b) 



where ks denotes the Boltzmann constant, and T is the temperature param- 
eter. 

In order to apply Eqs. (13), we require functions P{p, P) and p{P,p). From 
Eqs. (3), we find 



P(p,P) 



/M + m 



\M -m 
/M + m 

P{P:P)^( T7 J P 



p- 



/ 2M 

\M - m 
2m 



P 



,m - M J \m- M 
and, thus, Eqs. (13) take the explicit form 



(15a) 
(15b) 



Hp)^ 



M + m 



M -m 
M + m 



m- M 



dp0(p)$(P(p,P)) 
JdP^P)cf>{p{P,p)). 



(16a) 
(16b) 



As one can now easily verify by insertion, the Maxwell distributions (14) do in- 
deed satisfy Eqs. (16). Consequently, Eqs. (14) provide a pair of self-consistent 
stationary solutions for the non-relativistic collision kinematics, conditional on 
the information that a collision has occurred. 



3 Relativistic Brownian motions 



Since the integral equations (12) and (13) applied to the relativistic case as 
well, we can, in principle, proceed exactly analogous to the non-relativistic 
case. However, some purely technical difficulties arise due to the facts that (i) 
the relativistic collision kinematics is more complex than the non-relativistic 
one, and that (ii) the potential candidates for stationary distributions do not 
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allow for the solving Eqs. (12) and (13) analytically. Hence, after having speci- 
fied all required transformation formulae (Sec. 3.1), we will evaluate the PDFs 
on the left-hand-sides of the non-stationary integral equations (12) numerically 
^ , probing different types of candidate distributions (Sec. 3.2). 



3. 1 Kinematics of a single collision event 



In the relativistic case, the momentum and energy balance per (elastic) colli- 
sion can again be written in the form 



E + e^E + e, (17a) 
P+p^p + p, (17b) 

Compared with the non-relativistic case, the only difference is that we now 
use the relativistic expressions for momentum and kinetic energy, respectively. 
Specializing to the laboratory frame Sq (=rest frame of the heat bath), and 
using units such that the speed of light c = 1, we have [compare Eqs. (2)] 

P^MV^{V), p^mvyiv), (18a) 

where V and v denote the particles' velocities, and 



ry(v) ^ (l-V^) . (18c) 

Suppose we are given the information that a single collision has occurred. 
Then, solving Eqs. (17) for P, we find the explicit representations [compare 
Eqs. (3)] 



2?'or - (1 + t'o)P 
l-vl 



l-vi 

p(p,p)= ^''"';^\t'"^^ (19b) 



where the velocity 



P + P P + P . X 

-o=^ = T^=.o (19c) 



^ Since the integrals (12) are one-dimensional they can be solved numerically with 
e.g. the software package Mathematica [46]. 
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corresponds to the Lorentz boost from Sq to the center-of-mass frame (see 
App. A for details of the calculation). As one may easily check, in the non- 
relativistic limit case Eqs. (19) reduce to Eqs. (3). 

In order to be able to apply the integral criteria (12), we also need to determine 
the (partially) inverse transformations P{p,P) and p{P,p), respectively. In 
the non-relativistic case, this task was rather simple, see Eqs. (15). In the 
relativistic case, however, a bit of extra care is required. To illustrate this, let 
us first consider the momentum equation (19a) for Brownian particle. For any 
fixed value p with \p\ < oo, one finds 



pi^oo ^) = ^ ^20a) 

P.{p) :=^hm^P(p,P) = (20b) 

Hence, at finite p, the inverse transformation P{p, P) has a finite support, 
corresponding to the interval (without loss of generality, we will assume here 
and below that M > m) 



I{P) = [P-{P):P^ 

Hence, we find the following explicit form of the inverse transformation: 



P{p,P)^^^, Pel{p), pe (-00,00), (21a) 

with abbreviations 



Q = [m^ - + 2{p- P) p] [M^p - m^P + {p - P)] , (21b) 
R = 2{p- P) {m^ - M^) e E , (21c) 
5 = (m^ - M2)2 - 4(p - P) {M^p - m^P) . (21d) 

Note that the limits in Eqs. (20) correspond to the curves S — {). For later 
use, we also give the formal inversion of Eqs. (20): 



,^(Py.Jjll±ml±^^!l^JI^, (22a) 
^_^.^, J^^ + M^)P-JM^-n,)E ^ (22b, 

which allow us to rewrite Eq. (21) equivalently as 
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P{P,P) = ^^, Pe (-00,00), peb-(^),P+(^)]- (23) 

Equation (23) is in such a form that it can directly be inserted into the integral 
equation (12a). 

Finally, going through an analogous analysis for the heat bath particle yields 



p(P,p) = ^, pe [p_(P),p+(P)], Pe (-00,00), (24a) 



where 



? = [M2 -m^ + 2(P-p)P] [rri^P - M^p + rri^ {P - p)] , (24b) 
r^2{P-p){KP -m'^)eE, (24c) 
s = (m^ - M'^f -A{P-p) {m^P - M^p) , (24d) 



and 



' (P^ V ^(P ^ {m' + M')P + {M'-m')E 

p_(P) :=^hm^p(P,p) = ^ ^-^^ (25b) 

Note that the limits in Eq. (25) correspond to the curves s = 0, and their 
formal inversions can be defined by 



2m^ 

p_^yJnf^M-)i,-(M--ra-)l^ (26b) 
which allows us to rewrite Eq. (21) equivalently as 



p(P,p) = ^, pe (-00,00), Pe[P_(p),P+(p)]. (27) 

Equation (27) is given in such a form that it can directly be inserted into the 
integral equation (12b). 



3.2 Testing the integral criterion 



Analogous to the non-relativistic case, we aim to exploit the integral criteria 
(12) in order to determine the invariant momentum distributions. To this end. 
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we proceed as follows: For functions P{p,P) and p{P,p), given in Eqs. (23) 
and (27), we evaluate numerically the integrals (12): 



P 

iP) = 

Jp 



p+ip) 

-{p) 
p+ip) 



dp 



dP 



dP 



dP 

dp 



P-(p) 



0(p) $(P(p,P)), 
$(P) 0(p(P,p)). 



(28a) 
(28b) 



for a number specified grid-points P and p, respectively. The boundaries of 
the integration are chosen in accordance with the support intervals of the 
functions P(p, P) and p(P,p); cf. Eqs. (23) and Eqs. (27), respectively. We 
use different pairs of initial PDFs ($, 0), and test if 



$(p) ^ $(p), 0(p) ^ m (29) 

hold simultaneously. If the answer is positive, we conclude that the candidate 
functions ($, 4>) fulfill the stationarity criteria (13) and, thus, are invariant 
solutions for the relativistic coUision process. 

We next discuss our choice of the initial PDFs. Guided by the results of [26,27], 
we consider 



where rj e [0, 1] is a free parameter, and e and E denote the relativistic kinetic 
energies, respectively. By choosing the symmetric candidate distributions (30), 
we automatically specialize to the rest frame of the heat bath. For a given set 
of parameters (m, M, T, rj) the normalization constants 3\f,j(m) and ?sf^(M) are 
determined by the conditions 



dP $(P), 1 = / dp0(p). (31) 



Let us briefly recall how different values of 77 may arise in the context of the 
Langevin description of relativistic Brownian motions, as developed in [26]. 
Specializing to the rest frame of the heat bath, one can derive the following 
Langevin equation for the stochastic motion of a relativistic Brownian particle 
[26]: 



12 



dp 
lit 



' Er,{P) 

M 



1/2 



(32) 



Here, u is the friction parameter and ^ ordinary Gaussian white noise with 
amphtude D, i.e., ^ is characterized by 



{e(t)) = 0, (e(t)e(0) = 2D5(t-t')- (33) 

On the rhs. of Eq. (32) the noise ^ couples multiphcativcly to (a function 
of) the momentum coordinate P. Hence, different discretization rules may 
yield different stationary momentum distributions [47]. It is convenient to 
parameterize discretization rules as follows 



E,{P) = E[i^P{t) + {l-^)P{t + dt)i r/e [0,1], (34) 

with E{P) denoting the rclativistic kinetic energy. In this notation, e.g., r/ = 
corresponds to the HK post-point discretization rule [38-41], rj = 1/2 to the 
Fisk-Stratonovich mid-point discretization [34-36] and = 1 to Ito's pre-point 
discretization [32] . Introducing the temperature via the Einstein relation 



(35) 



the candidate PDFs in Eqs. (30) represent the stationary distributions asso- 
ciated with different values of 77 [26] . In particular, only for 77 = (HK rule) 
the standard relativistic Jiittner-Maxwell distribution is recovered. 



We are now in the position to discuss the numerical results. Figures 1 (a-c) and 
1 (e-f) show functions $(-P) and 0(p), respectively, as obtained for different 
values of rj. In each diagram the solid lines correspond to the initial momentum 
distributions from Eqs. (30). The triangles indicate the distributions resulting 
after the collision, $(P) and 0(]3), obtained by numerically integrating Eqs. 
(12) at 50 different values of P and p, respectively. As one readily observes, for 
r) — 1, corresponding to diagrams (c) and (f), the initial distributions remain 
invariant in the course the elastic collision process. Hence, according to these 
results, the modified Jiittner functions (30) with 77 = 1 are the relativistic 
analogs of the the non-rclativistic Maxwell distribution. For completeness, we 
mention that we have tested the integral criteria over a wide range of tem- 
perature and mass parameters and always found that only the distributions 
with 77 = 1 are invariant (all numerical integrations were performed with the 
function NIntegrate of the computer algebra program Mathematica [46]). 
Furthermore, we note that the normalization was conserved with high accu- 
racy during the numerical integration, i.e., the numerically found distributions 
($,0) remained normalized to unity with high accuracy. 
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Fig. 1. Initial distributions (solid line) and numerical solutions (triangles) of 
Eqs. (28) for the momentum PDFs of the Brownian particle (a-c) and the heat 
bath particles (d-f). As one readily observes, see diagrams (c) and (f), only for 
ij = 1 the initial distributions are left invariant by the elastic collision process. 



4 Summciry and discussion 



We have studied a simple microscopic model for ID non-relativistic and rela- 
tivistic Brownian motions. It was our main objective to identify the (invariant) 
stationary momentum distributions for both Brownian and heat bath parti- 
cles on the basis of the underlying microscopic collision processes. To this end 
we have formulated two integral criteria, which relate the initial momentum 
distributions to the momentum distributions resulting after an elastic collision 
of a Brownian particle with a heat bath particle (Sec. 2.3). The assumptions 
(postulates) underlying derivation of the integral criteria can be summarized 
as follows: 

• validity of the standard kinematic conservation laws, 

• occurrence of the collision event, 

• independence of the initial momenta. 

It was then demonstrated that, under these assumptions, the integral criteria 
do correctly reproduce the Maxwellian distributions as the stationary solutions 
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for the non-relativistic Brownian motions (Sec. 2.4). 

Subsequently, the integral criteria were applied to the relativistic case (Sec. 3.2). 
Here, we found that the standard Juttner distributions, corresponding to rj — 
in Eqs. (30), are not stationary with respect to the integral criteria; i.e., given 
the information that the collision has indeed occurred, the standard Jiittner 
distributions do not remain invariant in the course of an elastic relativistic 
collision. Instead, the invariant distributions are given by modified Jiittner 
functions, corresponding to rj — 1 in Eqs. (30). This result is quite surpris- 
ing, since initially we had expected that the invariant solutions are given 
by standard Jiittner functions with 77 = 0. However, as known from earlier 
work [26], modified Jiittner distributions with rj ^ may also appear as sta- 
tionary solutions in the ID relativistic Langevin theory, with the value of 
77 depending on the discretization scheme that is used. In this context, we 
mention a recent paper by Lehmann [48], who argues that Jiittner's original 
approach [42] is non-covariant. Furthermore, we note that the invariant solu- 
tion with f] = 1 can also be interpreted as a simple exponential (canonical) 
distribution) with respect to the Lorentz-invariant volume element of momen- 
tum space, d^p/p^ = d^p/{p'^ + m'^Y^'^, where D is the number of spatial 
dimensions [49]. 

Due to the fact that our results are based on only three basic assumptions, 
there is very little freedom for modifications such that one could hope to re- 
cover the standard Jiittner distributions as invariant solutions. Hence, accord- 
ing to our opinion, this problem deserves further consideration in the future. 
For example, as the next step, it would be desirable to perform similar studies 
for simple 2D and 3D models ^ . Then the kinematics of a single collision pro- 
cess becomes more complicated, because - even for the simplest hard-sphere 
models - momentum may be redistributed in different directions (depend- 
ing on the impact parameter). Hence, if one wishes to formulate analogous 
integral criteria for identifying the stationary (invariant) 2D/3D momentum 
distributions, one will have to include additional equations taking into account 
the cross-sections ^ . However, if it should turn out that deviations from the 
standard Jiittner-Maxwell distribution persist in higher space dimensions as 
well, then this might be of relevance for calculating relativistic corrections 
in high-energy physics [50] and astrophysics (e.g., to the Sunyaev-Zeldovich 

^ Very recently, the authors were informed by L. O. Silva that Marti et al. have 
implemented a 3D relativistic Montc-Carlo model for charged particle collisions, 
and that their numerically obtained equilibrium solutions correspond to modified 
Jiittner distributions with rj = 1 as well (conference poster by M. Marti, R. A. Fon- 
seca, L. O. Silva "A colHsion module for OSIRIS", P5.015, 32nd European Physical 
Society Conference on Plasma Physics, Rome, June 2006). 

^ We note, the applicability of simple kinematic models as discussed here is, in prin- 
ciple, limited to situations where high energies quantum processes, as e.g. creation 
and annihilation of particles, can be neglected. 
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effect [51,52]). 

The authors are grateful to J. Chluba, A. Fingerle, and, in particular, to S. 
Hilbert for numerous helpful discussions. 
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A Relativistic collisions 



Our goal is to validate Eqs. (19). To this end, consider two particles having 
positions x and X and relativistic momenta p and P, as defined by Eq. (18), 
with respect to (wrt.) the ID inertial lab-frame Eq. Assume that the particles 
collide at, say, t — 0, which implies that: 

(i) if v{t) > V{t) at i < 0, then x{t) < X{t) or, alternatively, 

(ii) if v{t) < V{t) at t < 0, then x{t) > X{t). 

In order calculate the momenta after the collision, p and P, as functions of the 
initial momenta, p and P, it is convenient to perform a Lorentz transforma- 
tion to the center-of-mass frame. Before doing this, we briefly establish some 
notations. We define (l-l-l)-momenta wrt. Eq by 

p = (e,p), V-(E,P). (A.l) 

Assume that some inertial frame E' moves with velocity vq in Eq, then the 
Lorentz transformation matrix A^g can be parametrized as follows 




(A.2) 

Its inverse is obtained by replacing vq with —vq, i.e.. 



A;o^ = A-o = ^=^ r (A.3) 

For example, given *P = {E,p), the corresponding (1-1-1) -momentum vector 
wrt. E', denoted by — {E' , P'), is obtained via matrix multiplication 

^' = A.0 ^- (A.4) 
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Center- of- mass frame. — In the following let us assume that E' is an 
inertial center-of-mass frame of the colliding particles; i.e., in E' we have, by 
definition. 



p' + P' = 0. (A.5) 
This condition determines the Lorentz transformation parameter as 



^0 = — TT = = ^0- (A.6) 

Furthermore, for elastic collisions, the energy and momentum balance become 
particularly simple in E': 

E' = E', e' = e', (A. 7a) 

P' = -P' =p' = -p', (A.7b) 

where, as before, hat-symbols refer to the momenta after the collision. It is 
convenient to express Eqs. (A. 7) in matrix form, e.g., writing 



«P' = (7'q3', p'^a'p', (A.8a) 
where the momentum transfer matrix a' of the elastic collision is defined by 



(x 



-1 



= (7'-^ (A.8b) 



Lab frame. - It is now straightforward to convert the above results (A. 8) to 
the laboratory frame by applying the corresponding Lorentz transformations; 
e.g., for the Brownian particle we find 



^ = K'l ^' = A,-; a' ^' = A,-^ a' A,, ^ = a qj, (A.9a) 



where 



^ ^0 \^ 2vq -(l+^;o) 

is the momentum transfer matrix wrt. Eq, and vq is given by Eq. (A.6). Anal- 
ogously, we find for the heat bath particle 
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An explicit evaluation of Eqs. (A. 9a) and (A. 10) yields 



1 ({l + vl)E-2voP' 
2voE-{l + vl 





1-vi ( 












~l-vi ( 



which contains the desired result, cf. Eqs. (19). 

Finally, we also calculate the inverse momentum transfer matrix 



a-' = (A-^ a' = A"^ a'-' = A"^ a' A 



allowing us to write 











-T 


-vi ( 








[pI 


-T 





{l + v^o)e-2vop' 



(A.IO) 



(A.lla) 
(A.llb) 



(A.12) 



'o \2voe - {l + v^)p 



(A.13a) 
(A. 13b) 
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